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Abstract 

A Majorana neutrino is characterized by just one flavor diagonal electromagnetic form 
factor: the anapole moment, that in the static limit corresponds to the axial vector charge 
radius (r^). Experimental information on this quantity is scarce, especially in the case of 
the tau neutrino. We present a comprehensive analysis of the available data on the single 
photon production process e + e~ — > vv^f off Z-resonance, and we discuss the constraints 
that these measurements can set on {r\) for the r neutrino. We also derive limits for 
the Dirac case, when the presence of a vector charge radius (r v ) is allowed. Finally, we 
comment on additional experimental data on scattering from the NuTeV, E734, CCFR 
and CHARM-II collaborations, and estimate the limits implied for (r\) and (r v ) for the 
muon neutrino. 
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1 Introduction 



Experimental evidences for neutrino oscillations [1, 2, 3, 4] imply that neutrinos are the first 
elementary particles whose properties cannot be fully described within the Standard Model 
(SM). This hints to the possibility that other properties of these intriguing particles might 
substantially deviate from the predictions of the SM, and is presently motivating vigorous 
efforts, both on the theoretical and experimental sides, to understand more in depth the detailed 
properties of neutrinos and of their interactions. In particular, electromagnetic properties of 
the neutrinos can play important roles in a wide variety of domains such as cosmology [5] and 
astrophysics [6, 7], and can also provide a viable explanation to the observed depletion of the 
electron neutrino flux from the Sun [8, 9, 10, 11, 12, 13]. 

The electromagnetic interaction of Dirac neutrinos is described in terms of four form factors. 
The matrix element of the electromagnetic current between an initial neutrino state Vi with 
momentum pi and a final state Vj with momentum pj reads [14, 15] 

(v?(Pj)\J™\v?(Pi)) = ™ 3 T°(q 2 )u t (1) 
rJV) = (q\ - qj) [V D (q 2 ) - A D (q 2 ) l5 ] + ta^[M D (q 2 ) + E D (q 2 ) l5 ] 

where q = Pj—pi, and the (ij) indexes denoting the relevant elements of the form factor matrices 
have been left implicit. In the % — j diagonal case, M D and E D are called the magnetic and 
electric form factors, that in the limit q 2 = define respectively the neutrino magnetic moment 
fj, = M D (0) and the (CP violating) electric dipole moment e = E D (0). The reduced Dirac 
form factor V D (q 2 ) and the neutrino anapole form factor A D (q 2 ) do not couple the neutrinos 
to on-shell photons. For % — j and in the q 2 = limit they are related to the vector and axial 
vector charge radii (rj>) and (r\) through 1 

(r 2 v ) = -6V D (0); (r 2 A ) = -6 A D (0) . (2) 

In the following, even when q 2 ^ we will keep referring to the reduced Dirac form factor 
and to the anapole form factor as the vector and axial vector charge radius. A long standing 
controversy about the possibility of consistently defining gauge invariant, physical, and process 
independent vector and axial vector charge radii [16] has been recently settled [17, 18, 19, 20]. 
The controversy was related to the general problem of defining improved one-loop Born ampli- 
tudes in SU(2) x U(l) for four fermion processes, like for example e + e~ — > //. If one tries to 
take into account one-loop vertex corrections by defining improved effective couplings, one finds 
that gauge invariance cannot be preserved unless, together with other one-loop contributions, 
W box diagrams are also added to the amplitude. However, box diagrams connect initial state 
fermions to the final states, and thus depend on the specific process. Due to the absence of 

1 The vector charge radius is defined as the second moment of the spatial charge distributions (vy) = 
J r 2 pv(r)df where pv(r) is the Fourier transform of the full Dirac form factor q 2 V D (q 2 ). The axial vector 
charge radius can be defined in a completely similar way. 
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a neutrino-photon coupling at the tree-level, the problem is even more acute when trying to 
define the charge radius as a physical, process independent property, intrinsic to neutrinos. In 
[17] it was realized that for neutrino scattering off right handed polarized fermions, the W box 
diagrams are absent to begin with, and thus no ambiguity arises. This suggested a way to derive 
a unique decomposition of loop contributions that separately respects gauge invariance, and 
from which a process independent charge radius could be defined as an intrinsic property of the 
neutrino. Furthermore, in [18, 20] it was argued that the so-defined charge radius is a physical 
observable, namely its value could be extracted, at least in principle, from experiments. 

For Majorana neutrinos, in the non-diagonal case {y^ 1 ^ vf 1 ) and in the limit of CP 
invariance the electromagnetic interaction is described by just two form factors [14]. If the initial 
and final Majorana neutrinos involved in the process have the same CP parity, only Ejf(q 2 ) and 
Ajf (q 2 ) are non vanishing, while if the CP parity is opposite, the electromagnetic interaction is 
described by M^(q 2 ) and V^(q 2 ). Finally, in the diagonal Majorana case v^ 1 = vf 1 the only 
surviving form factor is the anapole moment A M (q 2 ). As discussed in [21], this last result can 
be inferred from the requirement that the two identical fermions final state in 7 — > v M v M be 
antisymmetric, and therefore it holds regardless of the assumption of CP invariance. 

In the SM the neutrino electromagnetic form factors have extremely small values [22] . Due 
to the left-handed nature of the weak interactions, the numerical value of the vector and axial 
vector charge radius coincide, and for the different v e , and v T flavors they fall within the 
range (j~v,a) ~ (1 — 4) x 10~ 33 cm 2 [17]. 2 However, since neutrinos do show properties that 
are not accounted for by the SM, it could well be that also their electromagnetic interactions 
deviate substantially from the SM expectations. 

In general, the strongest limits on the neutrino electromagnetic form factors come from 
astrophysical and cosmological considerations. For example the neutrino magnetic moments can 
be constrained from consideration of stellar energy losses through plasma photon decay 7 — > vv 
[23], from the non-observation of anomalous energy loss in the Supernova 1987A neutrino 
burst as would have resulted from the rapid emission of superweakly interacting right handed 
neutrinos [23], and from Big Bang nucleosynthesis arguments. In this last case, the agreement 
between the measurements of primordial Helium abundance and the standard nucleosynthesis 
calculations imply that for example spin flipping Dirac magnetic moment interactions should 
be weak enough not to populate right handed neutrinos degrees of freedom at the time when 
the neutron-to-proton ratio freezes out [5]. 

Since the charge radii do not couple neutrinos to on-shell photons, the corresponding in- 
teractions are not relevant for stellar evolution arguments. However, in the Dirac case, right 
handed neutrinos can still be produced through e.g. e + e~ — > vr&r, and therefore constraints 
from the Supernova 1987A as well as from nucleosynthesis do apply. They yield respectively 
I (r 2 ) I < 2 x 10- 33 cm 2 [24] and I (r 2 ) I < 7 x 10- 33 cm 2 [25]. 3 

2 These values are obtained in the q 2 = limit, and decrease with increasing energies with a logarithmic 
behavior. 

3 In the SM with right handed neutrinos the vr cannot be produced through the charge radius couplings, 
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However, if neutrinos are Majorana particles, they don't have light right-handed partners, 
and the previous constraints do not apply. In this case, in particular for the r neutrino, an 
anapole moment corresponding to an interaction even stronger than electroweak could be al- 
lowed. In the early Universe such an interaction could keep v T in thermal equilibrium long 
enough to experience a substantial reheating from e + e~ — > v T v T annihilation. We have inves- 
tigated to what extent this reheating could affect the Universe expansion rate and change the 
predictions for primordial Helium abundance. As we will discuss in section 2, we have found 
that even an interaction one order of magnitude stronger than electroweak would hardly affect 
Helium abundance at an observable level. 

We conclude that constraints on the Majorana neutrino axial charge radius can be obtained 
only from terrestrial experiments. The present laboratory limits for the electron neutrino are 
—5.5 x 1CT 32 < (r\(u e )) < 9.8 x 10~ 32 cm 2 [26] . 4 Of course in the Dirac case these limits 
apply to the sum (r v ) + (r A ) as well. Limits for the muon neutrino have been derived from 
scattering experiments [27, 28]. They are about one order of magnitude stronger than for the 
electron neutrinos, and will be discussed in section 4. Due to the fact that intense v T beams 
are not available in laboratories, to date no direct limits on {r\{y T )) have been reported by 
experimental collaborations. However, under the assumption that a significant fraction of the 
neutrinos from the sun converts into is T , by using the SNO and Super-Kamiokande observations 
the limit |(r^(i/ T ))| p$ 2 x lCT 31 cm 2 has been derived [29]. A limit on the v T vector charge radius 
(Dirac case) was derived by analyzing TRISTAN data on the single photon production process 
e + e~ — > uuj [30]. The same data can be used to constrain also the anapole moment for a 
Majorana u T , and therefore we have included TRISTAN measurements in our set of constraints. 

In the next section we will briefly analyze the possibility of deriving constraints on the 
Majorana neutrino axial charge radius from nucleosynthesis. In section 3 we will study the 
bounds on the tau neutrino charge radius implied by the TRISTAN and LEP experimental 
results. In section 4 we will discuss the constraints on the muon neutrino charge radius from 
the NuTeV, CHARM-II, CCFR and the BNL E734 experiments. They result in the following 
90% c.l. limits: 

- 8.2 x 10~ 32 cm 2 < (r 2 A (u T )} < 9.9 x 10~ 32 cm 2 , (3) 
-5.2 x 10~ 33 cm 2 < (r\(u^)) < 6.8 x 10~ 33 cm 2 . (4) 

For (r 2 A {y e )) we could not find new experimental results that would imply better constraints than 
the existing ones [26]. We just mention that the Bugey nuclear reactor data from the detector 
module closest to the neutrino source (15 meters) [31] should imply independent limits of the 
same order of magnitude than the existing ones. 

since the vector and axial vector contributions exactly cancel. Therefore, the quoted limits implicitly assume 
that, because of new physics contributions, one of the two form factors dominates and no cancellations occur. 

4 These limits are twice the values published in [26] since we are using a convention for (r V A ) that differs for 
a factor of 2. 



4 



2 Nucleosynthesis 



In this section we study the possible impact on the primordial Helium abundance Y, of an axial 
charge radius large enough to keep a Majorana u T in thermal contact with the plasma down to 
temperatures T < 1 MeV. In this case the neutrinos would get reheated by e + e~ annihilation, 
and this would affect the Universe expansion rate. To give an example, if one neutrino species 
is maintained in thermal equilibrium until e + e~ annihilation is completed (T <C m e ) this would 
affect the expansion as Av = I — (4/If) 4//3 ~ 0.74 additional neutrinos. 

The amount of Helium produced in the early Universe is determined by the value of the 
neutron to proton ratio n/p at the time when the ne + <-> pv and nv <-> pe~ electroweak reactions 
freeze out. This occurs approximately at a temperature Tf Q pa 0.7 MeV [32, 33]. Apart for the 
effect of neutron decay, virtually all the surviving neutrons end up in A He nuclei. Assuming 
no anomalous contributions to the electron neutrino reactions, the freeze out temperature can 
only be affected by changes in the universe expansion rate, that is controlled by the number of 
relativistic degrees of freedom and by their temperature. If tau neutrinos have only standard 
interactions, at the time of the freeze out they are completely decoupled from the thermal 
plasma. However, an anomalous contribution to the process e + e~ <-> v T v T would allow the v T 
to share part of the entropy released in e + e~ annihilation. The maximum effect is achieved 
assuming that the new interaction is able to keep the v T thermalized down to Tf Q . The required 
strength of the new interaction can be estimated by equating the rate for an anomalously fast 
e + e _ <-> v T v T process Y Vt = (av)n e to the universe expansion rate IV = (87r p/3m 2 p ) 1 ^ 2 . In 
the previous formulas (av) is the thermally averaged cross section times relative velocity, n e pa 
0.365 T 3 is the number density of electrons, p ~ 1.66 g X J 2 (T 2 /mp) is the Universe energy density 
with g* pa 10.75 the number of relativistic degrees of freedom, and mp is the Plank mass. The 
thermally averaged cross section can be written as (av) ~ kG 2 v T 2 where G Vt pa (2ir 2 a/3) (r\) 
parametrizes the strength of the interaction and is assumed to be sensibly larger than the Fermi 
constant Gp, and k ~ 0.2 has been introduced to allow direct comparison with the SM rate 
(av) SM ~ 0.2G 2 F T 2 [32]. By setting V Vt = T v at T = T fo , we obtain G Vt w 13 x 10" 5 GeV- 2 . 
Therefore, to keep the v T thermalized until the ratio n/p freezes out, an interaction about ten 
times stronger than electroweak is needed. 

However, even in the presence of such a large interaction, Helium abundance would only be 
mildly affected. This is because at T pa 0.7 MeV e + e~ annihilation is still not very efficient, 
and the photon temperature is only slightly above the temperature of thermally decoupled 
neutrinos: {T^—T v )jT 1 pa 1.5% [32]. This induces a change in the primordial Helium abundance 
AV pa +0.04 (AT Ut /T„) which is below one part in one thousand. This effect could possibly 
be at the level of the present theoretical precision [34]; however, it is far below the present 
observational accuracy, for which the errors are of the order of one percent [35]. 
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3 Limits on v T vector and axial vector charge radii 



Limits on (r v ) and (r\) for v T can be set using experimental data on single photon production 
through the process e + e~ — > P1/7. In the following we will analyze the data from TRISTAN and 
the off-resonance data from LEP. These data have been collected over a large energy range, 
from 58 GeV up to 207 GeV. Given that form factors run with the energy, we will present 
separate results for the data collected below Z resonance (TRISTAN), between Z resonance 
and the threshold for W + W~ production (LEP-1.5), and finally for the data above W + W~ 
production (LEP-2). Due to the much larger statistics collected at high energy, a combined fit 
of all the data does not give any sizable improvement with respect to the LEP-2 limits, that 
therefore represent our stronger bounds. 

The SM cross section for the process e + e~ — > z/z/7 is given by [36] 

da ul/y 2a /n 



dxdy x(l - y 2 ) [\ 2/ 4 



i\ 2 x 2 y 2 
1-x + ' 



{N l/ a s (s',g v ,g A )+(J st (s')+a t (s')} (5) 



where a s corresponds to the lowest order s channel Z boson exchange with N v = 3 the number 
of neutrinos that couple to the Z boson. For later convenience in a s we have explicitly shown 
the dependence on the electron couplings gy = —1/2 + 2 sin 2 6w and g& = —1/2, where 6w 
is the weak mixing angle. The additional two terms a st and a t in (5) correspond respectively 
to Z-W interference and to t channel W boson exchange in v e production. The kinematic 
variables are the scaled photon momentum x = E^/E^m with E heilin = \fs/2, the reduced 
center of mass energy s' = s(l — x), and the cosine of the angle between the photon momentum 
and the incident beam direction y = cos 6* 7 . The expressions for the lowest order cross sections 
appearing in (5) read 

sGl \ (g 2 v + g\) Ml 
as[S) ~ 6vr (Mf - S ) 2 + Mfr| ' [ ) 

sG 2 F {gy + 9 A ) (Ml - s) M 2 
ast[) ~ 6tt (Ml- S ) 2 + Mlr| ' [) 

= #' (8) 

where Gf is the Fermi constant, a the fine structure constant, Mz and I 'z the mass and width 
of the Z boson. Few comments are in order. Eq. (5) was first derived in [36]. It holds at 
relatively low energies where W exchange in the t channel can be legitimately approximated as 
a contact interaction. This amounts to neglect the momentum transfer in the W propagator, 
and to drop the W-^y interaction, so that the photons are emitted only from the electron 
lines. While this approximation is sufficiently good at TRISTAN energies, to analyze the 
LEP data collected above Z resonance some improvements have to be introduced. We will 
use an improved approximation where finite distance effects are taken into account in the W 
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propagator, however we will still work in the limit of vanishing W-'y interactions. While strictly 
speaking the amplitude with photon attached only to the electron legs is not gauge invariant, 
the necessary contribution for completing the gauge invariant amplitude is of higher order in a 
leading log approximation [37], and for our scopes can be safely neglected. Finite distance W 
exchange effects can be taken into account in the previous expressions through the replacement 



v s t(s) 



where M w is the W boson mass, and 



<?st(s) ■ F st 



a t (s)-F t 




F st {z) = %l + z) 2 \og(l + z)-z(l + ^z) 



z u 
3 



Ft(z) = -2(l + z)\og(l + z) + z(2 + z) 



(9) 
(10) 

(11) 
(12) 



The contact interaction approximation is recovered in the limit z — * for which F sttt (z) — > 1. 

An anomalous interaction due to non-vanishing v T axial and axial vector charge radii can 
be directly included in (5) by redefining the Z boson exchange term in the following way: 



N v a s (s\ g v , g A ) -> (N v - 1) a s (s', g v , g A ) + cr s (s', g* v (s'),g A ). 



where 



9vW) = 
5 = 



gv - 

3G F 



1 - 



MI 



5, 



(r 2 v) + (r A ) 



(13) 

(14) 
(15) 



The substitution gv — > g v m (13) takes into account the new photon exchange diagram for 
production of left-handed v T . In the Dirac case, s-channel production of right handed v T through 
photon exchange must also be taken into account. This yields a new contribution that adds 
incoherently to the cross section, and that can be included by adding inside the brackets in (5) 
the term 



ctr(s') = 



5' 



b F /r/\2 



67T 

\f2nci 
3G F 



(ST , 

(r 2 y) ~ (r 2 A ) 



(16) 
(17) 



In the SM (r v ) = (r A ) and therefore there is no production of v R through these couplings. For 
a Majorana neutrino 5' — and (r v ) = 0, and thus the limits on anomalous contributions to 
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the process e + e" — > uuj translate into direct constraints on the axial charge radius (r\(v T )). 
Note that including anomalous contributions just for the v T is justified by the fact that for v e 
and the existing limits are generally stronger than what can be derived from the process 
under consideration. 

3.1 Limits from TRISTAN 

The three TRISTAN experiments AMY [38], TOPAZ [39] and VENUS [40] have searched for 
single photon production in e + e~ annihilation at a cm. energy of approximately y/s = 58 GeV. 
Anomalous contributions to the cross section for e + e~ — > z/z/7 would have been signaled by 
an excess of events in their measurements. Limits on the tau neutrino charge radius from the 
TRISTAN data have already been derived in [30]. In the present analysis, we include also 
the neutrino axial charge radius, and we give an alternative statistical treatment based on a 
X 2 -analysis and on the measured cross sections, rather than on the number of events observed 
combined with Poisson statistics as given in [30]. This puts the TRISTAN constraints on a 
comparable statistical basis with the LEP results discussed in the next section. 

TRISTAN data are collected in table 1. The number of single photon observed, including the 
SM backgrounds, was respectively 6 for AMY, 5 for TOPAZ, and 8 for VENUS. The numbers 
listed in the N ohs column in table 1 are the background subtracted events, that correspond 
to the measured cross sections a mcs given in the fourth column. We have found that our 
expressions for the cross section (5)-(12) tend to overestimate the Monte Carlo results quoted 
by the three collaborations. This might be due to additional specific experimental cuts besides 
the ones quoted in the last two columns in table 1 . In any case, the disagreements with 
the Monte Carlo results remain well below the experimental errors, and therefore we simply 
consider it as an additional theoretical uncertainty that we add in quadrature. In constructing 
the x 2 function, we use conservatively as experimental errors the upper figures of the three 
measurements. This is justified by the fact that the 7 — Z interference term arising from new 
physics is always sub dominant with respect to the square of the anomalous photon exchange 
diagram, and therefore new physics contributions would always increase the cross section. 

For a Majorana v T (5' = and (r 2 v ) = 0) the TRISTAN data imply the following 90 % c.l. 
limits: 

-3.7x 10~ 31 cm 2 < (r\(v T )) < 3.1 x 10" 31 cm 2 . (18) 

For the Dirac case, the associated production of right-handed states through <jr in (17) allows 
us to constrain independently the vector and axial vector charge radius. The 90 % c.l. are: 

-2.1 x 10~ 31 cm 2 < {r 2 V)A (u T )) < 1.8 x 10~ 31 cm 2 . (19) 

As we have already mentioned, strictly speaking the constraints just derived cannot be directly 
compared with the LEP constraints analyzed below, since the two experiments are proving 
neutrino form factors at different energy scales. Of course, since our limits are meaningful 
only to the extent that they are interpreted as constraints on physics beyond the SM, it is not 
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Table 1: Summary of the TRISTAN data: The center of mass energy and luminosity are 
given in the second and third column. The background subtracted experimental cross sections 
and the Monte Carlo expectations quoted by the three collaborations are given respectively in 
column four and five (in femtobarns), while the number of observed events after background 
subtraction is listed in column six. e is the efficiency of the cuts in percent units. The last 
two columns collect the kinematic cuts, with x = E 1 /E\, eaiU11 xt = xsin# 7 with d 1 the angle 
between the photon momentum and the beam direction, and y = cos# 7 . 





y/s [GcV] 




(j mcs [fb] 


a MC [fb] 






/ Eheam 


\y\ 






55 




34 




AA 


x>0.175 




AMY [38] 


57.8 


91 
56 
99 


for 
(x>0.125 

\y\ <o.7) 


34 
49 
49 


a q+3.7 a 
4 - z -2.6 


64 
58 
57 


x>0.175 
x>0.125 
x>0.125 


<0.7 


TOPAZ [39] 


58 


213 


07+58 
— 19 


54 


2.2±H b 


27.3 


x > 0.14 
x T > 0.12 


<0.8 


VENUS [40] 


58 


164.1 


ao n+ 45 - 3 

4Z - u -30.2 


36.4 


3.9±t;l c 


57 


x T >0.13 


<0.64 



a AMY observes 6 events in the 4 runs listed above (respectively 0, 2, 2, 2) with an estimated background of 
1.7 ± 0.3 events. The quoted value for iV bs has been derived from their background subtracted cross-section. 
6 TOPAZ observes 5 events, and expects 2.5^J;g from background. iVobs has been derived from their background 
subtracted cross section. 

c VENUS observes 8 events and expects 4.1^'y from background. They quote 3.9^2'g background subtracted 
vv^f events, which correspond to the cross section given in the fourth column. 

possible to make a sound guess of the form of the scaling of the form factors with the energy, 
which is determined by the details of the underlying new physics. However, if we assume a 
logarithmic reduction of the form factors with increasing energy as is the case in the SM, than 
we would expect a moderate reduction of about 0.65 when scaling from TRISTAN to LEP-1.5 
energies, and an additional reduction of about 0.75 from LEP-1.5 up to LEP-2 measurements 
at 200 GeV. 

3.2 Limits from LEP 

Limits on (ry) and {r\) can be derived from the observation of single photon production at LEP 
in a completely similar way. We stress that contrary to magnetic moment interactions that get 
enhanced at low energies with respect to electroweak interactions, the interaction corresponding 
to a charge radius scale with energy roughly in the same way than the electroweak interactions, 
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and therefore searches for possible effects at high energy are not in disadvantage with respect 
to low energy experiments. It is for this reason that LEP data above the Z resonance are able 
to set the best constraints on the vector and axial vector charge radius for the r neutrino. 

All LEP experiments have published high statistics data for the process e + e~ — > z/z/7 for cm. 
energies close to the Z-pole; however, due to the dominance of resonant Z boson exchange, these 
data are not useful to constrain anomalous neutrino couplings to s-channel off-shell photons. 
Therefore, in the following we will analyze LEP data on single photon production collected 
above Z resonance, in the energy range 130 GeV - 207 GeV. We divide the data into two 
sets: LEP-1.5 data collected below W + W~ production threshold are collected in table 2, while 
LEP-2 data, collected above W + W~ threshold and spanning the energy range 161 - 207 GeV 
are collected in table 3. 

3.2.1 LEP-1.5 

The ALEPH [41], DELPHI [42] and OPAL [43, 44, 45] collaborations have published data for 
single photon production at cm. energies of 130 GeV and 136 GeV. During the fall 1995 
runs ALEPH [41] and DELPHI [42] accumulated about 6 pb _1 of data for each experiment, 
observing respectively 40 and 37 events. In the same runs OPAL [43, 44] collected a little less 
than 5 pb _1 observing 53 events. In addition, OPAL published data also for the 1997 runs (at 
the same energies) [45] collecting an integrated luminosity of 5.7 pb _1 and observing 60 events. 

ALEPH reports two values for the cross sections at 130 GeV and 136 GeV, each based on 
2.9 pb _1 of statistics. They also quote the results of a Monte Carlo calculation of the SM cross 
section, that is in good agreement with the experimental numbers (and with our estimates). 
DELPHI combined together the statistics of both the 130 GeV and 136 GeV runs, however they 
present separate results for two different detector components: the High density Projection 
Chamber (HPC) covering large polar angles, and the Forward ElectroMagnetic Calorimeter 
(FEMC) covering small polar angles. Since DELPHI does not quote any Monte Carlo result we 
assign a bona fide 5% theoretical error for our cross section estimates. OPAL published two sets 
of data. The data recorded in the 1995 runs [43] were reanalyzed in [44], and correspond to 2.30 
pb" 1 collected at 130 GeV, and to 2.59 pb" 1 collected at 136 GeV. In the 1997 runs[45] 2.35 
pb _1 were collected at 130 GeV, and 3.37 pb _1 at 136 GeV. With a total integrated luminosity 
of about 28 pb -1 LEP-1.5 implies the following 90 % cl. limits: 

- 5.9 x 10" 31 cm 2 < (r 2 A (u T )) < 6.6 x 10~ 31 cm 2 (20) 

for the axial vector charge radius of a Majorana z/ r , and 

- 3.5 x 10" 31 cm 2 < (r^ A (u T )) < 3.7 x 10 -31 cm 2 (21) 

for the Dirac case. Let us note that, in spite of the much larger statistics, the limits from 
LEP-1.5 (20) and (21) are roughly a factor of two worse than the limits from TRISTAN in (18) 
and (19). The main reason for this is that at LEP-1.5 energies initial state radiation tends to 
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bring the effective cm. energy of the collision s' close to the Z resonance, thus enhancing Z 
exchange with respect to the new photon exchange diagram. 

Table 2: Summary of the ALEPH, DELPHI and OPAL data collected below W+W~ pro- 
duction threshold. ALEPH [41] and OPAL [44, 45] present separate results for two different 
energies, while DELPHI [42] combines together the data collected at 130 and 136 GeV. DEL- 
PHI presents separate data for two different detector components: the High density Projection 
Chamber (HPC) covering large polar angles, and the Forward Electromagnetic Calorimeter 
(FEMC) covering the forward regions. The kinematic cuts applied are given in columns eight 
and nine. Wherever a double error is listed, the first is statistical and the second is systematic. 



LEP-1.5 


y/s[GeW] 




a mcs [pb] 


a MC [pb] 


N obs 


e(%) 


E 1 [GeV] 


\y\ 


ALEPH 


130 


2.9 


9.6±2.0±0.3 


10.7±0.2 


23 


85 




















> io 


< 0.95 


[41] 


136 


2.9 


7.2±1.7±0.2 


9.1±0.2 


17 


85 






DELPHI 


















HPC [42] 


(133) 


5.83 


7.9±1.9±0.7 




20 


53* 


> 2 


< 0.70 


FEMC [42] 


(133) 


5.83 


6.0±1.9±0.6 




17 


43* 


> io 


0.83-0.98 


OPAL 


130 


2.30 


10.0±2.3±0.4 


13.48±0.22t 


19 


81.6 


x T >0.05 


< 0.82 


[44] 


136 


2.59 


16.3±2.8±0.7 


11.30±0.20 t 


34 


79.7 


or 

x T >0.1 


< 0.966 




130 


2.35 


11.6±2.5±0.4 


14.26±0.06 t 


21 


77.0 






[45] 








11.95±0.07 t 






x T >0.05 


< 0.966 


136 


3.37 


14.9±2.4±0.5 


39 


77.5 





* Estimated from the inferred experimental cross sections and measured numbers of events. 

t Calculated from the expected number of events as predicted by the KORALZ event generator. 



3.2.2 LEP-2 

Above the threshold for W + W~ production the four LEP experiments collected altogether 
about 1.6 nb _1 of data. The corresponding 24 data-points are collected in table (3). ALEPH 
[46, 47, 48] published data for ten different cm. energies, ranging from 161 GeV up to 209 
GeV. Data collected between 203.0 GeV and 205.5 GeV were combined together, they appear 
in the table as the 205 GeV entry, and the same was done for the data collected between 205.5 
GeV and 209.0 GeV that are quoted as the 207 GeV entry. DELPHI [49] published data col- 
lected at 183 GeV and 189 GeV, and gives separate results for the three major electromagnetic 
calorimeters, the HPC, the FEMC and the Small angle Tile Calorimeter (STIC) that covers 
the very forward regions, between 2° - 10° and 170° - 178°. In three papers [50, 51, 52] L3 
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reported the results obtained at 161 GeV, 172 GeV, 183 GeV and 189 GeV. While for most 
data points the agreement between our SM computation of the cross-sections and the Monte 
Carlo results is at the level of 5 % or better, we find that the L3 Monte Carlo results are up 
to 20% larger than our numbers, and this disagreement is encountered for all the four L3 data 
points. While we have not been able to track the reasons of this discrepancy, we have verified 
that the effects on our final results is negligible. OPAL published data for four different cm. 
energies [44, 45, 52]. For the data presented in [44, 45] we have estimated the Monte Carlo 
cross sections from the published numbers of events expected as predicted by the KORALZ 
event generator. The results agree well with our estimates. 
The 90 % c.l. limits implied by LEP-2 data read 

- 8.2 x 10~ 32 cm 2 < (r 2 A (v T )) < 9.9 x 10~ 32 cm 2 (22) 

for the Majorana case, and 

- 5.6 x 10~ 32 cm 2 < (r v>A (v T )) < 6.2 x 10~ 32 cm 2 (23) 

for a Dirac v T . 

These limits are about a factor of four stronger than the limits derived in [29] from the 
SNO and Super-Kamiokande observations and than the limits obtained in [30] from just the 
TRISTAN data. In Fig. 1 we depict the 90 % c.l. limits on (r v {v T )) and (r 2 A {v T )) for the Dirac 
case as derived from the LEP-2 data. The picture shows the absence of any strong correlation 
between (r v (v T )) and (r A (u T )). We stress that the possibility of bounding simultaneously the 
vector and axial vector charge radii stems from the fact that in e + e~ annihilation also the right- 
handed neutrinos can be produced, and they couple to the photon through a combination of 
(r 2 /) and (r A ) which is orthogonal to the one that couples the left-handed neutrinos. In contrast, 
neutrino scattering experiments do not involve the right handed neutrinos, and therefore can 
only bound the combination (r 2 ,) + (r\). 

Before concluding this section, we should mention that independent limits could also be 
derived from the DONUT experiment, through an analysis similar to the one presented in [54], 
and that yielded limits on the v T magnetic moment. We have estimated that the constraints 
from DONUT would be at least one order of magnitude worse than the limits obtained from 
LEP; however, it should be remarked that these limits would be inferred directly from the 
absence of anomalous interactions for a neutrino beam with an identified v T component [55]. 

4 Limits on vector and axial vector charge radius 

The NuTeV collaboration has recently published a value of sin 2 9\y measured from the ratio 
of neutral current to charged current in deep inelastic z/^-nucleon scattering [56] . Their result 

reads 

sin 2 9$ = 0.2277 ± 0.0013 ± 0.0009 (24) 
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Table 3: Summary of the ALEPH, DELPHI, L3 and OPAL experimental data, collected above 
W + W~ production threshold. The notation is the same than in (2). Wherever a double error 
is listed, the first is statistical and the second is systematic. 



LEP-2 


[GeV] 




a nics [pb] 


a MC [pb] 


N obs 


e(%) 


£ 7 [GeV] 


\y\ 


ALEPH 


161 


n.i 


5.3±0.8±0.2 


5.81±0.03 


41 


70 


x T > 0.075 


< 0.95 


[46] 


172 


10.6 


4.7±0.8±0.2 


4.85±0.04 


36 


72 


[47] 


183 


58.5 


4.32±0.31±0.13 


4.15±0.03 


195 


77 


x T > 0.075 


< 0.95 




189 


173.6 


3.43±0.16±0.06 


3.48±0.05 


484 










192 


28.9 


3.47±0.39±0.06 


3.23±0.05 


81 










196 


79.9 


3.03±0.22±0.06 


3.26±0.05 


197 








[481 

L W J 


200 


87.0 


3.23±0.21±0.06 


3.12±0.05 


231 


81.5 


x T > 0.075 


< 0.95 




202 


44.4 


2.99±0.29±0.05 


3.07±0.05 


110 










205 


79.5 


2.84±0.21±0.05 


2.93±0.05 


182 










207 


134.3 


2.67±0.16±0.05 


2.80±0.05 


292 








1 A I -T 1 T3 TXT 


















[49] 

L J 

HPC 


183 

i on 
lot) 


50.2 

1 r > 7 

154. ( 


1.85±0.25±0.15 

i onj_ n 1 tx_i_n 1 A 
l.o0±0.15±0.14 


2.04 
i.y i 


54 

1 A P. 

14o 


58» 


x > 0.06 


< 0.70 


FEMC 


183 
189 


49.2 
157.7 


2.33±0.31±0.18 
1.89±0.16±0.15 


2.08 
1.94 


65 
155 


54« 
50« 


x >0.2 
x <0.9 


>0.85 
<0.98 


STIC 


183 
189 


51. 4 
157.3 


1.27±0.25±0.11 
1.41±0.15±0.13 


1.50 
1.42 


32 
94 


f 

— f 

_ t 


x >0.3 
x <0.9 


>0.990 
<0.998 


L3 


161 


10.7 


6 75±0 91±0 18 

V7 • 1 / — 1 — V 7 • '.7 J 1 — V 7 • -L v7 


6.26±0.12 


57 


80.5 


> 10 
and 


< 0.73 




1 79 


i n 9 


a i O-L-n RQ-r-0 '\A 

U. IZiZlU.Oc/ZlU. ±^± 


n. «i _i_n i pi 


zLQ 


p,n 7 

ou. i 


H/T _ u 


fl Sfl-D Q7 


[51] 


183 


55.3 


5.36±0.39±0.10 


5.62±0.10 


195 


65.4 


> 5 


< 0.73 














and 




[52] 


189 


176.4 


5.25±0.22±0.07 


5.29±0.06 


572 


60.8 


£ T > 5 


0.81-0.97 


OPAL 


161 


9.89 


5.3±0.8±0.2 


6.49±0.08+ 


40 


75.2 


x T >0.05 


< 0.82 


[44] 


172 


10.28 


5.5±0.8±0.2 


5.53 ±0.08+ 


45 


77.9 


or 

x T >0.1 


< 0.966 


[45] 


183 


54.5 


4.71± 0.34±0.16 


4.98±0.02t 


191 


74.2 


x T >0.05 


< 0.966 


[53] 


189 


177.3 


4.35±0.17±0.09 


4.66±0.03 


643 


82.1 


x T >0.05 


< 0.966 



* The STIC Calorimeter efficiency varies between 74% and 27% over the angular region used in the analysis. 
" Estimated from the Monte Carlo cross sections and the expected numbers of events. 
^Calculated from the expected number of events as predicted by the KORALZ event generator. 
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-0.6-0.4-0.2 0.2 0.4 0.6 

(riK)) [10- 31 cm 2 ] 

Figure 1: Combined limits on {r v {v T )) and (r 2 A (v T )) for Dirac tau neutrinos derived from LEP-2 
data. The plot shows the Xmin + 2-71 contour, corresponding to 90 % c.l. 

where the first error is statistical and the second error is systematic. In order to derive limits 
on neutrino electromagnetic properties one should compare the results obtained in neutrino 
experiments to a value of sin 2 9w determined from experiments that do not involve neutri- 
nos. Currently, the most precise value of sin 2 6\y from non-neutrino experiments comes from 
measurements at the Z-pole and from direct measurements of the H^-mass [57]. In our nu- 
merical calculations we will use the value for sin 2 6 W obtained from a global fit to electroweak 
measurements without neutrino-nucleon scattering data, as reported in [56, 58]: 

sin 2 6 W = 0.2227 ±0.00037. (25) 

The effect of a non-vanishing charge radius can be taken into account through the replacement 
gv — > gv — 8 m the formulas for z/^-nucleon and z/^-electron scattering [59], where S is given 
in (15). Since there are no right-handed neutrinos involved, there is no effect proportional to 
5' and therefore only 5 oc {r v {v^)) + (r\{v^)) can be constrained. Upper and lower limits can 
be directly derived by comparing sin 2 0$ with the quoted value of sin 2 6 W from non- neutrino 
experiments. Since the results for neutrino experiments and the measurements at the Z-pole 
are not consistent at the la level, in the following equations (26)-(28) we will (conservatively) 
combine the errors by adding them linearly. 5 

From the NuTeV result (24) we obtain the 90 % c.l. upper limit: 

(r$W) + (^aM) < 7.1 x 10- 33 cm 2 . (26) 
5 Except for the CCFR data, which is consistent with the SM precision fits. 
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(b) 



(c) 



(d) 



Figure 2: 90 % c.l. limits on ((ry)+(r A )) for the muon neutrino derived from (a) E734 at BNL 
[28], (b) CHARM II [27], (c) CCFR experiment [60] and (d) from the NuTeV result [56]. 



A reanalysis of the E734 data on u^-e and v^-e scattering [28] yields the 90 % c.l. limits: 

- 5.7 x 10~ 32 cm 2 < (r v (^)) + (r A (v^)) < 1.1 x 10~ 32 cm 2 . (27) 

Note that in ref. [28] the E734 collaboration is quoting a lower limit about 3.6 times and an 
upper limit about 7.5 times tighter than the ones given in (27). This is because of various 
reasons: first of all, as was pointed out in [61], in [28] an inconsistent value for Gf was used 
that resulted in bounds stronger by approximately a factor of v^2. In addition, the errors were 
combined quadratically, which, due to the large negative trend in their data, resulted in a much 
stronger upper bound on {ry{v^) + {t\{v^)) than the one quoted here. Finally, our value of S 
is defined through the shift g v — > gv ~ $ of the SM vector coupling, consistently for example 
with the notation of [59], while the convention used by the E734 Collaboration [28] as well as 
by CHARM II [27] define 8 as a shift in sin 2 9w- This implies that our limits are larger for an 
additional factor of 2 with respect to the results published by these two collaborations. 
From the CHARM II neutrino-electron scattering data [27] we obtain at 90 % c.l.: 

- 0.52 x 10" 32 cm 2 < (r 2 ^)) + (r 2 A (^)) < 2.2 x 10" 32 cm 2 . (28) 

These limits differ from the numbers published by the CHARM II collaboration [27] not only 
because of the mentioned factor of 2 in the definition of 5, but also because the present value 
of sin 2 6\v [57] is smaller than the one used in 1995 in the CHARM II analysis. 
From the data published by the CCFR collaboration [60] one can deduce 

- 0.53 x 10~ 32 cm 2 < (r v {iy„)) + (r A (iy^)) < 0.68 x 10~ 32 cm 2 . (29) 
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The four limits discussed above are represented in fig. 2, that makes apparent the level of 
precision of the NuTeV result. By combining the upper limit from CCFR (29) and the lower 
limit from CHARM II (28) we finally obtain: 

- 5.2 x 1(T 33 cm 2 < (r v (is^) + (r 2 A (u^)) < 6.8 x 1(T 33 cm 2 . (30) 

It is well known that the NuTeV result shows a sizable deviation from the SM predictions 
[56], and as a consequence it also appears to be inconsistent (at the 90 % c.l.) with 5 = 0. 
In fact, strictly speaking their result {ry(u^)) + {r 2 A (u^)) = (4.20 ± 1.64) x 10~ 33 cm 2 (1 a 
error) could be interpreted as a measurement of (^y(^)) + ( r A( u n))- A vanishing value for 
( r v( l V)) + ( r A( u v)) becomes consistent with NuTeV data only at approximately 2.5 standard 
deviations. We should also mention that the fact that the NuTeV central value is very close to 
the range for (ry) at q 2 = quoted in the introduction should not mislead to think that a SM 
effect has been measured. In the SM the charge radius (r 2 (q 2 )) runs from its value at q 2 = 
approximately as \og(\q\ 2 / M^r) ■ In the NuTeV experiment the energy transfer is always > 20 
GeV [56], and therefore at the typical interaction energies of this experiment the value of the 
charge radius is expected to be smaller than its value in the static limit by at least a factor of 
ten. 

5 Conclusions 

This work stems from the observation that if neutrinos are Majorana particles their axial charge 
radius (r A ) , that is the only permitted flavor diagonal electromagnetic form factor, cannot be 
constrained through astrophysical or cosmological observations. In section 2 we have discussed 
in some detail how it is not possible to derive useful constraints from nucleosynthesis and 
from the measurements of primordial Helium abundance. We have concluded that in order to 
constrain (r A ) we can rely only on the analysis of the results of terrestrial experiments. 

In section 3 we have presented a comprehensive analysis of the available off Z-resonance 
data for the process e + e~ — > vv^. We have used these data to derive limits for the axial vector 
charge radius of the r neutrino, as well as combined limits on the vector and axial vector charge 
radius in the case of a Dirac v T . These limits are largely dominated by the high statistics LEP-2 
data collected above W + W~ production threshold. 

We have also analyzed the bounds that can be derived for the muon neutrino from an analysis 
of neutrino scattering experiments. We obtained the most stringent limits by combining the 
CCFR z/^-nucleon scattering and the CHARM II //^-electron scattering results. No new limits 
were obtained for the electron neutrino; however, new experiments dedicated to the detailed 
study of electron (anti) neutrino interactions with matter, as for example the MUNU experiment 
at the Bugey nuclear reactor [62], should be able to improve existing limits by about one order 
of magnitude. 
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